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AN ANALYTIC APPROXIMATION TO 
VERY HIGH SPECIFIC IMPULSE & SPECIFIC POWER 
INTERPLANETARY SPACE MISSION ANALYSIS 


Craig Hamilton Williams* 

A simple, analytic approximation is derived to calculate trip time and performance 
for propulsion systems of very high specific impulse (50,000 to 200,000 seconds) 
and very high specific power (10 to 1000 kW/kg) for human interplanetary space 
missions. The approach assumed field-free space, constant thrust/constant specific 
power, and near straight line (radial) trajectories between the planets. Closed form, 
one dimensional equations of motion for two-bum rendezvous and four-bum round 
trip missions are derived as a function of specific impulse, specific power, and 
propellant mass ratio. The equations are coupled to an optimizing parameter that 
maximizes performance and minimizes trip time. Data generated for hypothetical 
one-way and round trip human missions to Jupiter were found to be within 1 % and 
6% accuracy of integrated solutions respectively ; verifying that for these systems, 
credible analysis does not require computationally intensive numerical techniques. 

INTRODUCTION 

Through countless millennia, mankind has ventured to unknown worlds in the pursuit of wealth, 
knowledge, freedom, adventure, and international prestige. To date, human exploration of space has been 
primarily driven by the pursuit of international prestige. Because of the high cost and long trip times associated 
with human interplanetary exploration, such endeavors are expected to remain unattractive to any potential 
sponsor other than die governments of the most financially capable and technologically advanced nations. 

While it is difficult to imagine a low cost solution to human interplanetary travel, far-term technically feasible 
concepts do exist that could significantly reduce the expected long trip times. This could obviate larger and 
more complex spacecraft, higher operations costs, inordinately long return on investment times compared to 
commercial ventures, and the unappealing sacrifice of the spacefarer’s personal time. Conventional chemical, 
and nearer-term electric and nuclear fission propulsion technologies would require multi-year round-trip times 
to even the closest planets. Although long duration missions have been the norm for unmanned probes, they 
represent a real hardship for human travel. More advanced propulsion technologies, however, could provide 
trip times comparable to some commercial terrestrial operations — several weeks, perhaps even days. As D. 
Cole pointed out in his 1959 paper on minimum trip time travel to the planets, the evolution of terrestrial 
propulsion systems has been from the inexpensive “minimum energy” sailing vessels to the more costly 
“minimum time” aircraft '. Advanced space propulsion systems could catalyze future interplanetary travel 
predicated on the private, not public, sector market. High specific impulse(I^,)/high specific power (a) systems 
could accomplish this by traversing almost straight line trajectories between planets due to their expected 
ability to generate large accelerations. These propulsive technologies, unfortunately, are well beyond near- 
term feasibility and require solving some formidable technical challenges. 
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Because of the technological immaturity and perceived remoteness of these advanced propulsion 
concepts, most human interplanetary mission studies have been predicated on trip time and propulsion system 
performance data generated by analysis tools designed to model either 1^-limited systems (high thmst) or a- 
limited systems (low thrust). As a result, high 1,^/high a technologies are infrequently studied and therefore 
have few analytic tools designed to model them. Fortunately, due to the nature of their operation, simple 
analytical techniques can provide reasonably accurate estimates of their mission travel times and performance 
capabilities. These simple, first-order analytical methods can greatly assist in focussing de\ elopment of 
requisite precursor technologies. Very high thrust propulsion systems could produce accclerJions greater than 
the local acceleration due to solar gravity at Earth’s orbit (0.6 E-3 g). The normally thoughUof conics of 
minimum energy trajectories degenerate into straight line transfers at these acceleration levels, A “field-free 
space” approximation can then be invoked to greatly amplify the usually complex orbital mechanics. Gravity 
losses and optimum steering concerns can be neglected without introducing too much error, obviating the need 
for computationally intensive, numerically integrated solutions to support a proof of concept analysis. As a 
result, an analytic closed-form solution can be derived that calculates minimized trip time and maximized 
performance capability in a single equation (or at most two). 

An extensive literature search surfaced three published approaches that solve only select portions of 
this problem for a constant thrust engine; though it is also worth noting that a few papers have also been 
written on the simpler to analyze, constant acceleration (variable thrust) device '. During the 1950's, H. 
Preston-Thomas of the National Research Council of Canada derived several one dimensional velocity 
equations 2 . Found just prior to the publication of this paper, Preston-Thomas’ work contains trip time 
relations, but sparse derivations and few equations of motion made assessment of the approach difficult 
Wolfgang Moeckel of the NASA Lewis Research Center developed a simple analytic approach using a one 
dimensional “flat solar system” model to calculate trip time using the quotient of velocity and distance during 
the 1 960's \ High and low thrust systems were treated separately; the high Oirust relations based on vehicle 
AV calculated by the classic rocket equation, the low ihrust relations based on the time integral o” acceleraion 
squared equation (“J’ parameter), modified to optimize engine-on time. Moeckel’s approach had the advantage 
of extreme simplicity with reasonably good accuracy. However, the segregation ofhigi . and low thrust systems 
into separate equations did not illustrate the explicit dependence of trip time as an function of only distance, I^,, 
and e in advanced very high 1,,/very high a systems. 

Professor Dennis Shepherd of Cornell University developed an alternate approach in the 1 960's 4 . 
Based in part on some earlier work by Ernst Stuhlinger, Shepherd’s approach used a one dimensional, finite 
bum model (continuous thrusting with no coast phases). He derived a closed-form equation for distance 
traveled in fielc -free space by an accelerating rocket at constant thrust in terms of 1^,, a, and trip time. The 
classic rocket equation served as the starting point, later incorporating an optimizing parameter to maximize 
performance and minimis© trip time. One advantage of Shepherd’s finite bum approach was a more 
conservative representation of one dimensional motion. More important, however , was the explicit 
appearance of l v with a in a single equation optimized for payload and trip time, thus providing the analyst 
with a convenient way of evaluating the relative effect of distance, 1^ , and a on trip time for any mission 
distance. The main limitation was that the relation was valid for one-way flyby missions only, where the 
spacecraft would not rendezvous with the destination planet. 

The problem solved in this paper is 8n extension of Shepherd’s work. Two-bum “rendf zvous” and 
four-bum “round trip” mission equations of motion are derived with an optimization condition imposed. Each 
leg of the trip consists of an acceleration bum from Earth’s heliocentric position followed immediately by a 
deceleration bum (with no intervening coast period) which terminates at the heliocentric position of the 
destination planet (with zero radial velocity with respect to the sun). The transfers are assumed to originate and 
terminate outside the effective planetary gravity wells. These one dimensional straight line approximations 
along the heliocentric radius vector closely resemble numerically integrated solutions (Figure 1). 
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APPROACH 


Review and Reformulation of the Solution to the One-Burn Flyby Problem 


Shepherd’s 4 approach will be briefly 
summarized. The basic rocket equation fora 
constantly vehieje is: 
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The mass balance relation is given by: 
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Rearranging terms of the mass equation yields: 
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Substituting for the propellant mass ratio into figure 1 : Continuously Thrusting Trajectories 
the rocket equation yields: 
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From the definitions of jet power and specific power, a Structure to propellant mass ratio can be defined: 
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Shepherd introduced a parameter he called the characteristic velocity V c , which upon substitution illows for a 
useful expression for the structure to propellant mass ratio: 
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thus: 
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Substituting Eq. (9) into Eq. (4) yields a relation for payload ratio in terms of A V, exhaust and characteristic 
velocities: 
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Shepherd illustrated through a series of plots the implications of such an expression, primarily the existence of 
an optimal 1^ (in which the payload is maximized) for a particular mission A V. He demonstrated this by 
differentiating the above expression with respect to the parameter c/V c and setting it equal to zero. The result 
is the optimization parameter: 
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Shepherd then took the classical distance equation for a single bum, continuously accelerating rocket 
(neglecting gravity losses): 
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and solved for total distance traveled in time T (assuming zero initial velocity): 


s = ^[(l-X)ln(l-X)^] = cT 
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By substituting the optimization parameter Eq. ( 11 ) for the exponential ten® in the denominator of Eq. ( 1 3), 
we have: 
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At this point, Shepherd concluded with discussion of the implications of this acceleration-only distance 
equation. 

Whit is of primary interest is to rewrite Shepherd’s single bum flyby mission distance equation into a form 
where trip time is explicitly a function c. distance, 1^,, and a. Expressing c as gj,,, and recalling the definition 
of V c , Eq.(l 4) can be rewritten as a quadratic in T: 
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Solving for trip time via the quadratic formula: 
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Where only die positive radical term is meaningful. 


Shepherd demonstrated that c/V c has limits of »0.505<=c/V c <=l .0, where the lower limit corresponds to a 
zero payload ratio and the upj .r limit corresponds to a payload ratio of unity. These limits are a manifestation 
of imposing the optimization parameter on the rocket equation. Since the payload ratio is a function of c/V c 
and AV/c, and A V/c is a function of only c/V c , specifying c/V c is equivalent to specifying a payload ratio. 

Sir ce c/V e is solely a function of T, 1,^ and o , payload ratio can replace one of these three variables. For 
example, it was convenient to express the limiting a=f(I^,S) for the zero payload ratio case: 
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Substituting Eq. ( 1 8) into Eq (1 6) and selecting the meaningful root results in nn equation for the 1 uniting trip 
tunc 
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sp 


( 19 ) 


Thus the minimum trip time for a single bum flyby mission with negligible payload and optimal a varies 
linearly with distance and inversely with 1^, The optimal a can then be calculated from Eq. (18), 


Derivation of the Solution to the Two- Burn Rendezvous Problem 


Unlike the one-bum flyby mission, the rendezvous mission would include a deceleration bum to 
reduce the sun-relative radial velocity to zero by the time of arrival at the destination planet. The rendezvous 
mission also presupposed the utilization of in situ resources at the destination pli net to refuel the propulsion 
system, since the vehicle would carry only enough propellants for a one way trip. A modified rocket equation 
to model the deceleration bum is first needed. Shepherd’s time dependent form of the constant mass 
flow/constant thrust rocket equaUon for the acceleration bum ( 0 < t <= t,) was: 
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Therefore, a new deceleration bum (t, <= t < t, ) velocity function can be modeled similarly as: 
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V, is the velocity at the turn around point This point is where fee vehicle terminates acceleration, rotates 1 80 
degrees, and commences the deceleration bum. The logarithmic tom, therefore, is a negative or “braking” 

A V. The distance and acceleration expressions for the deceleration phase are easily derivable. To integrate the 
logarithmic term, the following substitution is convenient: 
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Completing the integration results in an equation for deceleration distance as a function of time 
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Taking the time derivatives of the velocity equations result in the acceleration mid deceleration equations as a 
function of time; 
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A representative plot of the position, velocity, and 
acceleration expressions for an example mission are 
illustrated in figure 2. 


Continuous A cceleration/Deceleration Mission 
Constant Thnm, lsp, & Power 
Example: Isp=50,000 sec and Aipha= 100 kW/kg 



Figure 2: Position, Velocity, & Acceleration 
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Thus, the stage propellant ratios are equal 


M n 

X. = - X 

1 M 


M 


Pi 


2 


r X 


( 30 ) 


Adding the two distance expressions (Eqs. (12) and (24), evaluated at their respective total bum times) and 
substituting Eqs. (28) and (30) results in an equation for total distance traveled; 
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The total distance traveled, S, is the straight line distance between the orbits of the two planets. It may be set to 
the difference between the radii (shortest trip time) or any other distance up to the opposition of the planets. 

Note that the deceleration phase propellant ratio with respect to the initial mass can be written: 
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A relationship can be found relating the two bum times by noting that the mass flow rates are equal and that 
the total mission time is merely the sum of the two bum phases: 
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After replacing t, in Eq. (31) with Eq. (33), and noting that t, is the total trip time T, the logarithmic terms 
cancel and the distance equation reduces to: 
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This is quite different from the acceleration-only distance Eq. (13), which is a logarithmic function of A The 
implication is that attempting to substitute the optimization parameter Eq. (11) (with «Ls exponenual function of 
A V/c) into the distance equation, ns was done in the ope -hum flyby case, will not result in a simple function of 
distance, 1^, and a (since there will not be a clean cancellation of the intermediate variable A V/c) The 
distance and optimization equations for the two-burn rendezvous case (as well as the four=hum round trip case 
as wc fhati see) arc therefore coupled by the dependent vanabk A and must be solved by iteration 

It was found to be easiest to have A serve as the iteration parameter. Remember that the A used above was 
defined ns a stage propellant ratio It must first be replaced with the equivalent expression for 'he total 
propellant ratio, which is implicit in the optimizing pninmeler equation 
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This quadratic in X r , r has a solution of : 
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For the cases of interest here, where ~0.505<= c/V,<=1.0, only the negative root term is meaningful. Using 
this expression to replace A^ with A^,., in Eq. (34) yields: 
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The optimizing parameter Eq. ( 1 1 ) can be expressed as a function of A^, 1^, and a by recalling the definition 
of A from the rocket equation and the definition of c/V c as a function of 1^, a, and T (Eq. (8». Substituting 
these relations, the optimizing parameter can be rewritten as : 


^ total _j 

( 38 ) 

Eqs. (37) and (38) are sufficient to solve the two-bum rendezvous problem. By choosing an initial value of 
A^ and substituting into Eq (38), an inifiul trip time can be found (From Shepherd, an initial (maximum) 
A^, value of 0.79681 2 is suggested. *) The trip time can then be substituted into Eq. (37) to find a new value 
of A,^. Repeated iteration using the average of the A Uial values as the new A (oUl was found to converge quite 
rapidly in most cases. 
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Once a converged solution is found, several dependent variables of interest can easily be calculated The first 
bum propellant ratio and c/V„ can be calculated from Eqs (36 ) and 18) respectively. The structure, payload, 
and initial A final thnist-to-weipht ratios are then given by 
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As in the single bum acceleration case, payload ratios can be specified and the limiting 1^,’s and a’s can be 
calculated. The value of the parameter c/V t for each corresponding value of payload ratio must be found first. 
This can be accomplished by solving the rocket equation (rewritten in terms of mass ratios) and the 
optimization parameter simultaneously after specifying the desired payload fraction. These equations are 
coupled, but were found to be easily solved by iteration. Shepherd had provided an approximate value of 
—■■0^05-for-frte-zero-paylead-caser-Altemate-eHd-moFe-e^act- values- were-desired-for analysts purposes—For ■©%•-■ 
and 25% payload fractions, the values of c/V c were calculated to be 0.504976295 and 0.736886943 
.. respectively. Tfre vftlue of A .^. r pan then be found by inserting the optimization parameter Eq...( 1 1 ) into the.. . . 

rewritten rocket equation Eq. (10) (replacing the exponential term), substituting Eq. (36) for the equivalent 
AV/c term, and solving for The result is: 
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One* X is known, the optimum travel time cart bo cnsily calculated from Eq. (34 ) The optimum re can llx n 
be calculated from Eq. (8 ) At this point, the Uvo-bum problem and the parameters of interest arc solved 
These relation -i can be manipulated by band or easily assembled into a short computer routine 

Derivation of the Solution to the Four-Bum Round Trip Problem 

A similar net of equations were derived for a lbur-bum round trip mission, where the spue* -raft would cam nil 
of its propellant without the use of in situ propellant refueling, The equations were it mnd t< » be sitmlui In tli* • 
nvu-burn rendezvous mission Tbetwo-b'im rendezvou : mission can lie thought of ns tin: return leg of u 
round trip mission. Tliu velocity equations are of the same form as the two=bum mission, as an the 
acceleration equations Mince Uie outgoing and reluming, distances will be assumed to he equal 
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The relationships between trip times and X's must be found. Since the muss flow rates are constant: 
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Equating Eqs. (45) and (46), X, and k } are thus related by: 
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An expression for X's in terms of X ^ can be found since 


X 


total 






M, M t M, M, 


A,+A/l -X,)+X 3 (l -X,) 2 +A 3 (l -X 3 )(l -A,) 2 


(48) 


11 



Substituting Eq, (47) into Eq. (48) results in a quadratic expression for A, in terms of A tortl : 
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Only the negative radical term is meaningful. Notice that this first bum propellant fraction is half the 
expression found in the two-bum rendezvous case (Eq. (36)). 

Finally, the outgoing trip lime can be expressed as a fraction of the total round trip time by again noting that 
mass flow rate is constant: 
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Substituting Eq. (50) into the round trip distance Eq. (44): 
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As was done for the rendezvous case, this can be solved for A^ using Eq. (49). For convenience, the distance 
can be written in terms of S , rather than S ^ : 
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Using this equation and the optimization parameter Eq. (38), the round trip mission trip time can be calculated 
in the same way as was done for the rendezvous case. The c/V c , structure, and payload ratios are calculated the 
same as in the rendezvous mission (Eqs. (8), (39), and (40)). The optimum fixed payload (Eq. (43)) is reduced 
by a factor of one half as per Eq (49). The thrust-to-weight ratios are: 
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RESULTS 


Iterations on the trip time/optimization parameter equations were well behaved and converged 
rapidly. Although sample hand calculations were not found to be overly burdensome, a simple computer 
routine was written to calculate almost all of the following data. As an example, rendezvous and round trip 
r-iissions to Jupiter were calculated. After this data was generated, the same mission-spaee parameters were 

and establish the limits of validity of this approach. The position of Jupiter was assumed to be at minimum 
distance from Earth (4.203 AU) at departure. Round trip missions assumed no dwell time at the destination 
and no significant relative motion of the planets. Ranges of 1^, from SO, 000 to 200,000 seconds and a from 
1 00 to 1 000 kW/kg were representative of what a propulsion system must operate at to perform multi-week/ 
multi-month transfers to this nearest of major planets. Overall propulsion system efficiency was set to unity. 
Figure 3 illustrates rendezvous mission trip time as a function of both and a. A rendezvous trip time to 
Jupiter of six weeks, for example, will require propulsion technologies with I^'s of at least 50,000 seconds and 
a’s of at least 100 kW/kg. For constant 1^ , increasing a always decreases trip time since the engine jet power- 
out increases (and/or the si i ucture mass decreases). For constant a, increasing 1^, always lengthens trip time 
as the propulsion system oegins to act more like a low thrust electric device. 1^,’s beyond the 50,000 to 
1 00,000 seconds levels are r-xoessive for interplanetary AV’s mless a corresponding (significant) 
improvement in a also occurs. A propulsion system engineer, planning to conduct proof of concept 
experiments, can use these observations to readily trade-off work on one parameter for the other depending on 
the judged engineering difficulty. Propulsion technology research can thus be guided by the knowledge of the 
relative merit of pursuing greater thrust vs. greater power. 


Using the fixed payload mass fraction equations, payload contours can be superimposed onto the 1^, 
& o plots as a means of measuring performance. Two performance levels were calculated, a 0% (fastest trip 
time) and a 25% payload mass fraction. Other desired payload mass fractions between these two values can be 
approximated by interpolation. Greater values should not be estimated in this way due to the non-linearity of 
the function as c/V c increases. For trip times of approximately six '«eeks, the maximum payload mass fractions 
for I^’s between 50,000 ana 100,000 seconds were 3% to 25% respectively. To halve trip times (three 
weeks), 1^,’s would have to be doubled and a’s would have to be improved by an order of magnitude. 

Figure 4 illustrates the results for round trip missions, using similar ranges of 1^, and a. In the round 
trip mission, the spacecraft departing from Earth carries all the fuel required to perform the entire mission. By 
observation, the data illustrates that (like terrestrial travel) refueling ft a destination is a compelling way to 
operate. For e- -ample, a zero payload/ 50,000 second 1^, vehicle that refuels at a destination would take half 
the time (2 X 40 days) to perform the same mission as a same initial weight vehicle carrying all of its 
propellant from the start. This enhances the attractiveness of travel to the major outer planets due to then- 
abundant supply of potential fuels and propellants such as hydrogen, deuterium, and helium-3. These 
materials are plentiful in both free molecular form in planetary atmospheres and in bound forms on the surfaces 
of their moons. Major planets would not only be destinations, but transportation nodes as well, supplying 
resources for fuels/propellants and human consumption. 


Figure 5 and Table 1 contain additional rendezvous mission data: mass fractions, trip times, and 
thrust-to-weight data for 1^ = 50,000 seconds plotted as a function of a. For fixed 1^, as a increases, the 
amount of propellant needed increases at the expense of payload fraction. The structure ratio reaches a 
maximum at approximately 26% as Shepherd’s work had predicted 4 . As trip time is reduced by increasing a, 
payload ratio becomes vanishingly small. Thnist-to-weight ratios depart from near constant values as a 
increases. Figure 6 and Table 2 contain similar data for round trip missions, plotted as a function of 1^. 
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Round Drip Time (days) Rendezvous (one way) Trip rime (days) 


Example: Mission to Jupiter (Rendezvous) 

Trip Time vs. Isp & Alpha 
(Maximized payload ratio, minimized travel time) 



Example: Mission to Jupiter (Round Trip) 
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Figure 5. Selected Mission Data 

Example: Mission to Jupiter (Rendezvous) 


Mass Ratio , Trip Time, & Thrust -to-Weight vs. Specific Power 

Isp = 50,000 Sec 


Ismxrrurtl 




UproptUial — IpropsUm 


Tri^^m Piym^framcn 




Table 1. Selected Mission Data 


Specific 

Power 

(kW/kg) 

Time 

(days) 

Payload 

Ratio 

Isp 

Structure 

Ratio 

= 50,000 Sec 

Total 1st Bum 

Prop Prop 

Ratio Ratio 

c/Vc 

FAVt 

Initial 

Accel(g’s) 

FAVt 

Final 

Deceit's) 

10 

187.9 

0.528 

0.200 

0.271 

0.146 

0.860 

0.835E-3 

0.115E-2 

25 

97.7 

C 282 

0.260 

0.458 

0764 

0.754 

0771E-2 

0.500E-2 

50 

64.4 

0.128 

0.263 

0.609 

0.375 

0.657 

0.548E-2 

0.140E-1 

75 

52.0 

0.063 

0.246 

0.691 

0.444 

0.597 

0.769E-2 

0749E-1 

100 

45.3 

0.028 

0.228 

0.744 

0.494 

0.554 

0.950E-2 

0.371E-1 




Alpha = 100 kW/kg 




Specific 

Time 

Payload 

Structure 

Total 

1st Bum 

c/Vc 

FAVt 

FAVt 

Impulse 

(days) 

Ratio 

Ratio 

Prop 

Prop 


Initial 

Final 

(Sec) 




Ratio 

Ratio 


Ac'.el(R’s) 

Deceit's) 

50,000 

45.3 

0.028 

07.28 

0.744 

0.494 

0.554 

0.950E-2 

0.37 IE-1 

100,000 

79.3 

0.468 

0.219 

0.313 

0.171 

0.837 

0.457E-2 

0.665E-2 

150,000 

143.1 

0.758 

0.113 

0.129 

0.067 

0.934 

0.157E-2 

0.180E-2 

200,000 

236.5 

| 0.882 

0.057 

0.061 

0.031 

0.969 

0.595E-3 

0.634E-3 
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Payload Mass Radk® 










Figi re 6. Selected Mission Data 

Example: Mission to Jupiter (Round Trip) 


Mass Ratio , Trip Time, & Thrust -to-Weight vs. Specific Impulse 


Alpha = 25 kW/kg 




50,000 100,000 150,000 200.000 

Specific Impulse (Sec) 



F/WQntttfF/WtFmal 



Specific Impulse (Sec) 


TbpTtae 


Trip Tone 





Table 2. Selected Mission Data 


Specific Impulse (Sec) 


Isp = 200,000 Sec 


Specific 

Power 

(kW/kg) 

Time 

(days) 

Payload 

Ratio 

Structure 

Ratio 

Total 

Prop 

Ratio 

1st Bum 
Prop 
Ratio 

c/Vc 

F/Wt 

Initial 

Accel(g's) 

F/Wt 

Final 

Deceit's) 

25 

946. 1 

0.882 

0.057 

0.061 

0.015 

0.969 

0.149E-3 

0.158E-3 

50 

499.5 

0.788 

0.100 

0.112 

0.029 

0.943 

0.520E-3 

0.585E-3 

75 

349.7 

0.711 

0.132 

0.156 

0.041 

0.920 

0.103E-2 

0.123E-2 

100 

274.2 

0.647 

0.158 

0.195 

0.051 

0.900 

0.165E-2 

0204E-2 

250 

134.8 

0.407 

0236 

0.357 

0.099 

0.812 

0.614E-2 

0.956E-2 

500 

84.8 

0228 

0266 

0.507 

0.149 

0.724 

0.138E-1 

0280E-1 

750 

66.7 

0.140 

0264 

0.595 

0.182 

0.666 

0.207E-1 

0.510E-1 

1000 

57.1 

0.090 

0255 

0.655 

0206 

0.624 

0266E-1 

0.771 E-l 

Specific 

Impulse 

(Sec) 

Time 

(days) 

Payload 

Ratio 

Alpha = 25 kW/kg 

Structure Total 1st Bum 

Ratio Prop Prop 

Ratio Ratio 

c/Vc 

F/Wt 

Initial 

Accel(g’s) 

F/Wt 

Final 

Decel(g's) 

50,000 

181.1 

0.028 

0228 

0.744 

0.247 

0.554 

0238E-2 

0.927E-2 

100,000 

317.2 

0.468 

0219 

0.313 

0.086 

0.837 

0.114E-2 

0.166E-2 

150,000 

572.6 

0.758 

0.113 

0.129 

0.033 

0.934 

0.392E-3 

0.450E-3 

200,000 

946.1 


0.057 

0.061 

0.015 


0.149E-3 

0.158E-3 
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Table 3. Equation Summary 



1 Bum 
Flyby 

2 Bum 
Rendezvous 

4 Bum 
Round! up 

Total 

Distance 

(13) 

(r 

(52) 

Optimized 
Trip Tune 

(16) 

(37)&(38) 

(52)&<38) 


Position 

Velocity 

Acceleration 

Acceleration 

Phase 

(12) 

(20) 

(25) 

Deceleration 

Phase 

(24) 

(21) 

(26) 


To check the validity of these results, a series of integrated computer runs were performed using/an 
existing, high fidelity interplanetary trajectory optimization routine 5 A. few of the results of these runs art'i 
superimposed onto the data generated and are shown in the last p! Figures 5 and 6. It was found that it 

these high power levels and with accelerations well in excess of the . ocal acceleration due to solar gravityitrip 
time results agreed within 1% for rendezvous and 6% for round trip missions. The larger error in round nip 
times was noticeable in cases where travel times where a considerable fraction (or greater) of one Earth year. 

In these cases, the integration program significantly altered the outgoing trip to accommodate the change in 
relative position of the Earth upon return. Since the approach outlined in this paper does not account for 
change in file relative positions of the planets implicitly, closer agreement for round trip missions is limited to 
the time scales similar to those of rendezvous missions. Table 3 provides an overall summary of pertinent trip 
time relations and equations of motion for the analyst 

CONCLUSIONS 

A simple, analytic approximation was shown to provide a means to readily calculate trip time and 
performance of propulsion systems of very high 1^ & a for human interplanetary space missions. Simultaneous 
solution of two equations was shown to be all that was necessary without having to resort to complex computer 
programs to integrate trajectories. The simplifying assumptions of one dimensional motion, field-free space, 
constant thmst/I tp /a, and near straight line (radial) trajectories between the planets permitted data of sufficient 
accuracy to be generated, providing insight into the relationships between distance, 1^ , and o. An optimizing 
parameter was included to maximize performance and minimize trip time. Example data was provided for 
rendezvous and round trip missions to Jupiter, illustrating trip times and payload mass fractions for a wide 
range of 1^, and a. Results indicated that a propulsion system technology must be capable of 1^,’ s of at least 

50.000 seconds and a’s of at least 100 kW/kg to travel to Jupiter in 6 weeks (rendezvous) and still allow for a 
3% payload mass fraction. Overall results were shown to agree within 1% and 6 % (rendezvous and round trip, 
respectively) of data independently generated, verifying that for preliminary analysis, a technology planner need 
not resort to large, high fidelity computer programs that require expert operators in order to accurately 
characterize the propulsion system parameter-space. Comparison of data between two-bum rendezvous and 
four-bum round trip missions illustrated the extreme attractiveness of in situ refueling at the destination planet. 
Few propulsion technologies have been identified that a»e expected to be capable of providing such 
demanding I^’s and a’s. Inertial confinement fusion systems may be capable of I^’s between 50,000 to 

270.000 seconds and a’s as large as 100 kW/kg *. Antiproton-catalyzed fusion and antiproton annihilation 
are other potential space propulsion concepts that are expected to operate at or beyond this regime. 
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NOTATION 


A acceleration (ft/sec**2) 

a specific power (ft lbf7(sec lbm)) 

c exhaust velocity (ft/sec) = gj^ 

F thrust (Ibf) 

g gravitational acceleration: 32.1739 (ft/sec**2) 

g* conversion constant: 32.1 739 (lbm/lbf)(ft/sec**2) 

I, specific impulse (lbf sec/lbm) = thrust per mass flow rate 
X propellant mass fraction (My/Mj) 

M mass (lbm) 

T( overall propulsion system efficiency 

P jet power (ft lbfrsec) 

S distance (ft) 

t time (sec) 

T trip time (sec) 

V 0 characteristic velocity (ft/sec) 

V velocity (ft/sec) 

AV velocity increment (ft/sec) 

W weight (lbf) 

Subscripts 

0,1,23,4 propulsion system burns & phases of flight 
f final 

i initial 

p propellant 

pay useful payload 

s structure (including tankage & power supply) 
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